First Passage Properties of Molecular Spiders 
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Molecular spiders are synthetic catalytic DNA-based nanoscale walkers. We study the mean first 
passage time for abstract models of spiders moving on a finite two-dimensional lattice with various 
boundary conditions, and compare it with the mean first passage time of spiders moving on a one- 
dimensional track. We evaluate by how much the slowdown on newly visited sites, owing to catalysis, 
can improve the mean first passage time of spiders and show that in one dimension, when both ends 
of the track are an absorbing boundary, the performance gain is lower than in two dimensions, when 
the absorbing boundary is a circle; this persists even when the absorbing boundary is a single site. 



I. INTRODUCTION 

Natural molecular motors play an important role in bi- 
ological processes that are critical for the functioning of 
living organisms; they are the source of most forms of mo- 
tion in living beings [l|-|3| . In addition to naturally occur- 
ring molecular motors, several synthetic molecular mo- 
tors have been designed |4T-fl2| . Our work is inspired by a 
particular type of synthetic molecular motors — molecular 
spiders. 

Molecular spiders [Tj, [lj] are synthetic nanoscale walk- 
ers which consist of a rigid, inert chemical body to 
which multiple flexible legs are attached. The legs are 
deoxyribozymes — enzymatic sequences of single-stranded 
DNA that can bind to and cleave complementary strands 
of a DNA substrate. When many such substrates are at- 
tached to a surface, a leg can move between substrates, 
cleaving them and leaving behind product DNA strands. 
Products can be revisited by a leg, but they cannot be 
cleaved again (Fig. [T]). The leg cleaves, and then de- 
taches, more slowly from a substrate than it detaches 
from a product. The number of legs, and their lengths, 
can be varied, and this defines how a spider moves on the 
surface, i.e., its gait. 

Mathematical models of molecular spiders at various 
levels of abstraction have been proposed and studied. 
Antal and collaborators introduced the first abstract 
model of molecular spiders, and studied the motion of 
a single spider on a one-dimensional track. They investi- 
gated the movement of spiders with various numbers of 
legs and various gaits over products only [151 ] and showed 
that such spiders are equivalent to a regular diffusion; the 
diffusion constants were computed for some gaits. Subse- 
quently they introduced substrates and took into account 
that cleavage and detachment from substrates together 
take more time than the detachment from products [la ] , 
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FIG. 1: (Color online) Molecular spider moves over a surface 
of single-stranded DNA substrates. It has several flexible de- 
oxyribozyme legs. When a leg detaches, it cleaves a substrate 
strand, turning it into a shorter product DNA strand. The 
leg can re-attach, but the bond will be weaker. 



showing that this difference in residence time and the 
presence of multiple legs, together, bias a spider's mo- 
tion towards fresh substrates when it is on a boundary 
between substrates and products. This important proj 
erty was also observed experimentally 



141. In Ref. h 



we showed that spiders move supcrdiffusively for long 
periods of time. Samii et al. investigated various gaits 
and numbers of legs [18|, UM,, emphasizing the possibil- 
ity of detachment from the track. In Refs. [2(J, [21| we 
studied the behavior of multiple spiders continuously re- 
leased onto a ID track. In a model with more physical de- 
tail [22[ , we showed that spiders can move against a force 
applied to the body. Models of spiders in two dimen- 
sions have also been studied. In Ref. [23| we investigated 
how fast several spiders with various gaits can locate a 
small number of targets placed on a small fixed-size two- 
dimensional lattice. In Ref. [2J| Antal and Krapivsky 
evaluated the diffusion constant and the amplitude de- 
scribing the asymptotic behavior of the number of visited 
sites for a single spider with various gaits placed on an 



infinite square lattice. Analytical results regarding the 
asymptotic behaviors (limit theorems, transience, recur- 
rence, and rate of escape) of spiders have been derived in 
Refs. [Hii]. In Refs. [27|,[28| the behavior of spiders in 
random environments was studied. Rank et al. showed 
that several spiders, each placed on a separate ID track 
and connected to a single cargo particle move it faster 
and remain superdiffusive longer than a single spider on 
a single ID track [H]. 

Here we study first passage properties of an abstract 
theoretical model of molecular spiders. Our model is a 
direct extension of the model introduced in Refs. (la Ha]. 
Although it is inspired by real molecular spiders, the 
model can also be applied to a wider class of random 
walkers that exhibit properties similar to spiders. Par- 
ticularly, we investigate how the various boundary condi- 
tions affect the spider's mean first passage time (MFPT) 
when it moves over finite one- and two-dimensional sur- 
faces. 

First passage properties of regular random walkers 
moving over discrete surfaces with various reflecting 
and absorbing boundaries have been extensively stud- 
ied [30, |31( . Here we show that the difference between the 
time the legs spend on visited and on unvisited sites re- 
duces the MFPT of two-legged spiders in various surface 
settings, and increases the MFPT of one-legged spiders 
(which behave like regular random walkers) . 

We start with the model of a two-legged spider 
that moves over a one-dimensional track with absorb- 
ing boundaries at both ends of the track. We found that 
for this surface the cleavage rate significantly affects the 
MFPT, and for any track length there exists an optimal 
cleavage rate. Next, we study an extension of the ID 
model to 2D, i.e., the mean first passage time of a two- 
legged spider to a circle, where the spider starts in the 
center. For this surface we determined that the cleav- 
age rate gives the spider an even greater advantage over 
a regular random walker. The advantage persists even 
when the target is a single site, and thus is much harder 
to find. In this second 2D model the circle is a reflecting 
boundary, its center is an absorbing boundary, and the 
spider starts from various distances from the center. 



II. MODEL DETAILS 

Our model is a modification of the model we used in 
Ref. [23|. It also can been seen as a direct extension of 
the AK model of the spiders on a plane [24j that takes 
various boundary conditions into account. 



A. Motion and Surface 

In our model a single spider moves over finite one- or 
two-dimensional regular lattices. A spider has k legs. It 
moves by detaching a leg from its site on the lattice and 
reattaching it to a new site. Only one leg can detach at 



any given time, so a spider cannot detach all of its legs 
to leave the lattice. Each site can be occupied only by 
one leg at a time. There is a restriction on the max- 
imum distance between any two legs S (the gait), and 
each leg can move to one of the nearest neighboring sites 
(2 sites in ID and 4 sites in 2D; a diagonal step is not 
allowed) with equal probability as long as the move does 
not violate one of the constraints above. Here we use 
only two types of spiders. First, a spider with k = 1; 
this spider is equivalent to a regular random walker. The 
parameter S does not affect this spider since it has only 
one leg. Second, a spider with k = 2 and 5 = 2; this 
spider is exactly the bipedal Euclidean spider with max- 
imal separation 2 [2ij . When such a spider is placed on 
a one-dimensional lattice, the model becomes equivalent 
to that of Ref. [l7J . 

Two types of sites can be present on the lattice, sub- 
strate and product. A leg detaches at rate 1 from a 
product, and at rate r from a substrate, where normally 
r < 1. Reattachment is instantaneous. When a leg leaves 
a substrate, that substrate is transformed into a product. 
Initially all sites are substrates; so any unvisited site is 
always a substrate, and a visited site is always a prod- 
uct. For r < 1 the legs act differently when they are 
on visited versus unvisited sites, as described above. On 
the other hand, when r = 1, legs act effectively the same 
whether they are on the products or substrates, and it 
becomes irrelevant if a site is visited (a product) or un- 
visited (a substrate). Thus, spiders with r < 1 can be 
seen as having memory, since they react differently to 
visited and unvisited sites, and spiders with r = 1 can be 
seen as having no memory, since they do not make this 
distinction. 



B. Boundaries and Starting Positions 

We study the first passage time of spiders moving 
over one-and two-dimensional regular lattices with var- 
ious boundary conditions and initial configurations. In 
each case we are interested when the spider reaches the 
absorbing boundary. 

The boundaries effectively make all our surfaces finite. 
In one dimension we use a ID track of length 2x. The 
spider starts its movement from the middle of the track 
(the origin) , and absorbing boundary sites are located at 
each end of the track, i.e., each one is x sites away from 
the origin. 

In two dimensions the surface is bounded by a circle of 
radius x. Here we study three types of boundaries and 
initial spider positions. First, the circle is an absorbing 
boundary, and the spider starts from the center of the 
circle. Second, the circle is a reflective boundary, while 
the center is an absorbing boundary, and the spider can 
start from any site on the circle. Third, we study the 
case when radius x is fixed and the spider starts y sites 
away from the target site at the center. 

In all these settings x effectively defines the size of the 



surface. And in all those settings, except the last, we 
study the dependence of MFPT on x, i.e., (t(x)). In the 
last case we study the dependence of MFPT on y, i.e., 

(r(y))- 



III. SPIDER ON A ID TRACK 

Background: Transient Superdiffusivity of 
Spiders 



Meaningfulness of MFPT in The Studied 
Settings 



A process to determine if MFPT is a valid character- 
istic of the first passage behavior was given in Ref. (32| ■ 
Similarly to that, we assess the meaningfulness of the 
mean first passage time in all studied settings. For every 
setting we estimate the distribution P(uj) of the random 
variable u = ti/(ti+T2), where t± and t<± are first passage 
times of two independent spiders. Values of w close to 1/2 
indicate that spiders act similarly in a particular setting. 
When u is close to or 1, the process is not uniform, and 
MFPT is not a good measure of actual behavior. Distri- 
bution P(oj) can have three distinct shapes: unimodal 
bell-shaped, bimodal M-shaped, and plateau-like, almost 
uniform behavior. Bell-shaped form with a maximum at 
uj = 1/2 indicates that MFPT can be considered as a 
valid measure of the first passage times of individual spi- 
ders. M-shaped form with two peaks close to and 1, 
and local minimum at 1/2 indicates that MFPT is not a 
good measure of the first passage time of individual spi- 
ders. The plateau-like shape with zero second derivative 
at co = 1/2 separates the two above cases. Just as in 
Ref. [32|], to quantify the shape of P(u>) we fit P(lu) to 
the model x^ 2 + c\uj + C2 for 0.05 < uj < 0.95. The sign 
of x indicates the shape of P(co). In cases when % < 0, 
the distribution is bell-shaped; x > shows that the dis- 
tribution is bimodal, M-shaped; and x = indicates that 
the distribution is almost uniform. 



D. Simulation 



Combining the states of the spider and the surface 
gives us a continuous-time Markov process for our model. 
We use the Kinetic Monte Carlo method [33[ to simulate 
many trajectories of the Markov process for every in- 
stance of the parameter set. In each case we record the 
first passage time to the absorbing boundary. 



For Section lrnl we simulated 210 different r rates for up 
to a distance of 10000 using 21504 traces of the Markov 
process. For Section llVl we simulated 100 different r rates 
for up to a distance of 1000 using 20000 traces. For 
Sections IV Al we simulated 18 different r rates for up to 
a distance of 250 using 20000 traces. For Section IV Bl we 
simulated 15 different r rates for up to a distance of 100 
using 20000 traces. 



Antal and Krapivsky analytically obtain the mean time 
(T(n)) to visit n sites in one dimension [15|, [Tg|. Eq. Q] 
shows how (T(n)) depends on the value of r for a two- 
legged spider (k = 2). 



(T(n)) 



31 



22 



(1) 



The parameter r affects the leading asymptotic behav- 
ior of (T(n)). Eq.HJ also from Ref. 0, gives (T(n)) for 
a one-legged spider 



(T(n)) = 



n{n — 1) 



2r 



(2) 



Eq. [2] demonstrates that for the one-legged spider 
(k = 1), the leading term of (T(n)) does not depend 
on the parameter r. The r only slows the one-legged spi- 
der down by increasing the sub- leading term. Eq. [5] and 
Eq. [I] also show that in the absence of memory (r = 1) 
the one-legged spider is faster than the two-legged spi- 
der. Thus these two properties, which separately make 
the walkers slower, surprisingly improve the performance 
of the spiders when combined together. This happens 
because difference in residence time between visited and 
unvisited sites biases multi-legged spiders towards un- 
visited sites when they find themseves on the boundary 
between visited and unvisited sites. 

Using Kinetic Monte Carlo [33| simulations of the 
Markov process we showed the unanticipated result that 
spiders of the AK model with k = 2, S = 2, and r < 1 
move superdiffusively over a significant span of time and 
distance before eventually slowing down to move diffu- 
sively [17] . 

Superdiffusive motion can be described using mean 
square displacement of a walker as a function of time. 
The mean squared displacement is given by Eq. [3l where 
d is the number of dimensions, and D is an amplitude. 



(X 2 ) = 2dDt° 



'a = stationary 

< a < 1 subdiffusive 
a = 1 diffusive 

1 < a < 2 superdiffusive 

a = 2 ballistic or linear 



(3) 



From Eq. [3] we derive the condition for the walker to 
be superdiffusive at time t. The walker is moving instan- 
taneously superdiffusively [34| at a given time t if 



a(t) 



d(ln((X 2 )(t))) 
d(\n(t)) 



> 1. 



(4) 



In [l7(, we showed that each spider process goes 
through three different phases of motion defined by its 
value of a. Initially spiders arc at the origin, and must 
wait for both legs to cleave a substrate before they start 
moving at all. So when t < \jr the process is essen- 
tially stationary (the initial phase). After the spiders 
with r < 1 take several steps, they show a sustained pe- 
riod of supcrdiffusivc motion over many decades in time 
(the superdiffusive phase). Finally, as time goes to in- 
finity, all spiders will approach ordinary diffusion with 
awl (the diffusive phase); the spiders mainly move 
over regions of previously visited sites, which makes the 
value of r less relevant. 



B. Mean First Passage Time 

We measure the mean first passage time, (r(x)), where 
x is the absolute distance of the walker from the origin 
on a one-dimensional track. At that point the walker 
is absorbed by the boundary. Fig. [2] shows the initial 
configuration of the track where the walker is positioned 
in the middle, and the absorbing boundaries arc shown 
as stars. The length of the track is 2x. 
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FIG. 2: Initial configuration of the track. All sites are initially 
substrates. Absorbing boundary is represented by stars. 

According to the results of our numerical simulations 
the MFPT of both one- and two-legged spiders is pro- 
portional to x 2 . Eq. [5] shows the leading and sub- leading 
terms of (r(x)). 



(t(x)) ~ A\(r)x 2 + ai(f)x. 



(5) 



The amplitude A\ of the leading term describes the 
asymptotic behavior of the MFPT. For the one-legged 
spider we found that A\ does not depend on r. For the 
two-legged spider A\ increases with r, and approaches 2 
when r = 1. As r approaches zero A± ps 1.4. Fig.[3]shows 
the amplitude A\ for one- and two-legged spiders for 182 
r values less than 1. 

For both one and two-legged spiders the amplitude a\ 
of the sub-leading term decreases monotonically and ap- 
proaches in the absence of memory (r = 1). The ampli- 
tudes A\ and a\ show that the parameter r only increases 
the MFPT of one-legged spiders. For two-legged spiders 
varying r can decrease their MFPT. 

In Section IIII Al we recalled that when r < 1, the 
AK-model walkers go through three different regimes 
of motion — the initial, superdiffusive, and the diffusive 
stage. For lower r values the initial slow period is longer 
than for higher r values, but subsequently the superdif- 
fusive period is longer and faster. For travel over shorter 
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FIG. 3: Amplitude Ai from Eq.[S]as a function of the cleavage 
and detachment rate r. The data are for one-legged (k — 1) 
and two-legged (k = 2) spiders on a one-dimensional track. 



distances the initial period is more important and thus 
larger r values result in lower first passage times. For 
travel over longer distances the superdiffusive period is 
important and smaller r values give better results. Thus 
for every particular distance there is an optimal value of 
r that minimizes the MFPT. For example, for distance 
2000 the spider with r = 0.05 is faster than the other 
(sampled) r values; but for distances 4000 and longer the 
spider with r = 0.01 is faster. 

We estimated r opt (x) for tracks of various lengths (up 
to 10000) through simulations of 210 r values. The re- 
sults are shown in Fig. 2) The figure also shows the (t(x)) 
that corresponds to r opt (x) for each distance, i.e., the 
minimum (t(x)) achievable by varying the parameter r. 
The optimum r monotonically decreases with distance. 
Smaller r values create a stronger bias towards unvisited 
sites, but they make individual steps slower when the 
spider is on the boundary. Fig. @] shows that for better 
performance on shorter distances faster steps are more 
important than the bias towards unvisited sites, whereas 
for longer distances the bias dominates the MFPT. 

We also estimate r opt (a;) by fitting. First, we assume 
that A\{r) and a\{r) have the functional form of Eq. [51 
and find the constants c\ to cs by fitting Eq. [5] to the 
estimates of A\ and a±. 



Ai(r) = (cir + c 2 )/(c 3 r + c 4 ) 
ai(r) = {c 5 r + c 6 )/(c 7 r + c 8 ) 



(6) 



Then, we substitute the results into Eq. [5] and find 
when (t(x)) is minimized by extracting the derivative 
of (t(x)) with respect to r. The predicted r opt (x) is also 
drawn in Fig.UJ The derivation also shows that r op t(a;) ~ 
x^ 1 / 2 (see the Appendix for details). 

Our estimates of the indicator \ f° r various r and x 
values show that in the setting of the one-dimensional 
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FIG. 4: Optimal r values (r op t(x)) for various distances, and 
the corresponding {t(x)). Optimal r values are obtained in 
two ways; first, by fitting the data for (t(x)) into the model of 
Eq. [5] and amplitudes A\ and a± into the model of Eq. [6] and 
second, by simulating 210 values of r and choosing those that 
correspond to the lowest (t(x)). The estimation of r op t(x) 



shows that r opt (a;) 



-1/2 



track x is always negative, independent of the values of 
the parameters x and r. Thus the MFPT is a meaning- 
ful, reliable measure of the first passage time of spiders 
moving on a one-dimensional track. 



IV. SPIDER ON A 2D PLANE WITH A 
CIRCULAR ABSORBING BOUNDARY 

A direct extension of a one-dimensional track with 
length 2x to two dimensions is a set of sites bounded 
by a circle of radius x. Fig. [5] shows the initial config- 
uration of the surface with the spider positioned in the 
middle. 
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FIG. 5: Initial configuration of the 2D surface. All sites are 
initially substrates. The spider starts at the center of the 
circle. The absorbing boundary is shown as a circle; as soon 
as either leg crosses the circle the target is considered to be 
found, and the experiment stops. 

Using numerical simulations we found that, similarly 
to the spiders in one dimension, the MFPT of both one- 



and two-legged spiders is proportional to x 2 . However, 
the sub-leading terms arc much closer to the leading 
terms, and therefore are more important for estimating 
the MFPT. This implies that in two dimensions spiders 
(especially those with very small cleavage and detach- 
ment rate, r < 0.05) approach the asymptotic behavior 
especially slowly. Eq.[7]shows the leading and sub-leading 
terms of (t(x)). 



-(*)} 



A 2 {r)x 2 + a 2 (r)x 2 /lnt 
A 2 (r)x 2 +a 2 (r)x 2 /(\ntf- S8 



(7) 



The correlation (lni) -088 for the one-legged spider is 
unusual and slower than for the two-legged spider. Some- 
what similar effects were observed in Ref. [24| | in the es- 
timation of the mean squared displacement. As in one 
dimension, the amplitude A 2 of the leading term does 
not depend on r for the one-legged spider and is propor- 
tional to r for the two-legged spider. The sub-leading 
term's amplitude a 2 monotonically decreases with r in 
both cases. Fig. [6] shows the amplitude A 2 for the one- 
and two-legged spiders for 60 r values that are less than 
1. 
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FIG. 6: Amplitude A2 from Eq.[7]as a function of the cleavage 
and detachment rate r. The data is shown for the one-legged 
(k = 1) and two-legged (k — 2) spiders moving over two- 
dimensional surface with a circular absorbing boundary. 

For one-legged spiders, as in one dimension, lower r val- 
ues only increase the MFPT. The comparison of Figs. [S] 
and [3] shows that for two-legged spiders A 2 is affected 
more strongly by r than A\ . The A 2 of the two-legged 
spider in two dimensions even intersects the A 2 of the 
one-legged spider. The A 2 starts at 2 when r = I and 
decreases towards ps 0.5 as r approaches zero. 

Similar to the one-dimensional case, spiders with r < 1 
start more slowly than the no- memory spider with r = I ; 
then they move faster, and finally they slow down and 
approach regular diffusion. But in contrast to spiders on 
a ID track, the r values that correspond to fastest times 



are higher, and the MFPT of spiders with very small r 
values (less than 0.1) approaches the MFPT of spiders 
with r = 1 very slowly. However, the transition towards 
the diffusive stage happens more slowly compared with 
ID. 

We estimated r opt (x) for circles of various sizes (up to 
1000) through simulations of 100 r values, choosing the 
fastest ones for each distance. Comparison of the results 
is shown in Fig. [7] The figure also shows the fastest (r(x)) 
that corresponds to r opt (x) for each distance. Similarly 
to the ID case, we also estimate r Q pt(x) by fitting. 
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FIG. 7: Optimal r values (r op t(x)) for various radii. Optimal 
r values are obtained in two ways; first, by fitting the data 
for (t(x)) into the model of Eq. [7] and amplitudes Ai and ai 
into the model of Eq. [6j and second, by simulating 100 values 
of r and choosing those that correspond to the lowest (r(x)). 

Our estimates of the indicator \ for various r and x 
values show that in the setting of the two-dimensional 
circle with absorbing boundary at the perimeter x is al- 
ways negative, and this does not change with the values 
of the parameters x and r. Just as in ID, the MFPT 
is a meaningful measure of the first passage time of the 
individual spiders moving on a two-dimensional lattice. 



study how (r(x)) is affected by the parameter r, and how 
it is affected by x in (1) and y in (2). 



A. 2D Circle of Variable Radius With Target in 
the Middle and Spider Starting from the Boundary 

Fig. [5] shows the initial configuration of the surface 
where the spider is positioned on the contour and the sin- 
gle target site (absorbing boundary) is shown as a star. 
Since the reflecting boundary is a circle and the absorb- 
ing boundary is a single site in the center, all starting 
positions of the same distance from the target site are 
equivalent. In our simulations we pick a fixed position 
on the contour as a starting position for the spider. This 
position is shown as two small solid circles in Fig. [8] 
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FIG. 8: Initial configuration of the surface. All sites are ini- 
tially substrates. The absorbing boundary is shown as a star. 
The spider starts at the periphery. The circle is a reflecting 
boundary and thus steps outside of the circle are not allowed. 

We measure the first passage time of the walker from 
the contour of the surface to its center for various radii 
x. 

Interestingly, the MFPT asymptotically grows slightly 
faster than x 2 . The sub-leading term also grows faster 
than in the circular absorbing boundary case. Eq. [5] 
shows the leading and sub- leading terms of (t(x)). 



V. SPIDER ON A 2D PLANE WITH A 

CIRCULAR REFLECTING BOUNDARY AND AN 

ABSORBING BOUNDARY IN THE CENTER 

When we change the contour of the 2D surface to be 
a circular reflecting boundary instead of an absorbing 
one, and place a single target site in the center, spiders 
with memory (r < 1) still have an advantage over spiders 
without memory (r = 1). We consider two cases: (1) 
when the radius x of the circle is variable, and spider 
starts from any point on the contour, and (2) when the 
radius x is fixed, and the distance x between the starting 
position and the target is variable. Since both cases are 
circularly symmetric, all starting positions with the same 
distance from the center arc equivalent. In both cases we 



(r(x)) ~ B 2 (r)x 2 + b 2 (r)x 



1.5 



(8) 



Fig. [9] shows the amplitude B 2 for the two-legged spi- 
ders for 20 different r values between 0.1 and 1.0. 

Values of r smaller than 1 give the spider an advan- 
tage, similar to the 2D configuration discussed above, 
but now even for shorter radii x. The typical time-line 
of this process can be characterized as follows. First, 
the spider starts moving, it eventually visits many sites 
without finding the target site, and leaves behind many 
smaller regions of substrates. Next, the spider starts to 
move only over visited sites more often, and eventually 
encounters regions of substrates of various shapes and 
sizes. At this stage, the spider with r = 1 will not be 
affected by the substrate regions, and will move just as if 
they were visited sites. On the other hand, spiders with 




FIG. 9: Amplitude B-z from the Eq. Has a function of the 
cleavage and detachment rate r. The data are shown for two- 
legged (k = 2) spiders moving over a two-dimensional surface 
with a circular reflecting boundary. The single site in the 
middle of the circle is an absorbing boundary. 



r < 1 will become biased to stay on the substrates and 
explore those regions more thoroughly; this will increase 
their chances of finding the target site, since it must be 
in one of those unvisited areas. This scenario can explain 
how spiders with r < 1 gain an advantage over spiders 
with r = 1 on this surface. 

The indicator x-, m this setting, grows with the param- 
eter x, and for x gi 17 x becomes close to 0. As in the 
previously described settings, x does not depend on the 
parameter r. Fig. [10] shows the dependence of x on x 
for r = 0.5. For small surfaces (small values of x), the 



measure of the first passage time of individual spiders. 
However, for larger x, close to values of x indicate that 
the shape of the distribution P(lS) is close to uniform, 
and thus the MFPT is not as good a measure of indi- 
vidual behavior as it is in the settings when the spider 
searches for the contour. It also shows that the possible 
paths that the spider can take to locate a single target 
are more diverse than paths that lead to the contour. 



B. 2D Circle of Fixed Radius With Target in the 
Middle and Spider Starting From Various Distances 

Fig. [TT] shows the initial configuration of the surface 
where the walker is positioned at a distance y from the 
target site, and the radius x is set to 100 and remains 
constant. 
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FIG. 11: Initial configuration of the surface. All sites are 
initially substrates. The absorbing boundary is shown as a 
star. The spider starts y sites away from the center. A circle 
of constant radius x = 100 is a reflecting boundary and thus 
steps outside of the circle are not allowed. 
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FIG. 10: Indicator \ shows how well the MFPT can describe 
the first passage times of the individual spiders to a single tar- 
get. The surface is two-dimensional with a circular reflecting 
boundary of radius x, and spider starts from the contour. 

X is negative, and thus in those cases MFPT is a good 



We measure the first passage time of the walker from 
the contour of the surface to its center for various dis- 
tances y. Fig. H21 shows a plot of the MFPT. 

The shape of the curves is asymptotically logarithmic; 
this shows that the initial position of the spider docs 
not significantly affect (r(y)) when x is fixed. Even if 
the spider starts closer to the target, there are too many 
possible paths to the target in 2D, of which one is ran- 
domly chosen. Many of them are very long and initially 
lead the spider far away from the target. 

The indicator x^ m this setting, is positive for y Js 80 
and decreases with the parameter y. For y ^ 70 as the 
surface configuration approaches the configuration dis- 
cussed in the Section TV A[ and x becomes negative, but 
as in the case of the Section IV Al it still remains close to 
0. Fig. IT51 shows the dependence of x on y for r = 0.5. 
The positive values of x f° r t ne smaller y values indicate 
that MFPT is not a good measure of the first passage 
times of the individual spiders, and paths that lead the 
spider to a target are very diverse and vary significantly 
in their lengths. 
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FIG. 12: Mean first passage time to a single point in the center 
of the circle with a fixed radius. The boundary is reflecting, 
and spider starts at various distances from the target. 
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FIG. 13: Indicator \ shows how well the MFPT can describe 
the first passage times of the individual spiders to a single 
target. The surface is two-dimensional with a circular reflect- 
ing boundary of a fixed radius x = 100, and spider starts at 
the distance y from the target site. 



VI. COMPARISON OF SPIDER 
PERFORMANCE IN THE STUDIED SETTINGS 

It is interesting to compare the advantage spiders with 
r < 1 (i.e., with memory) enjoy over those with r = 
1 (i.e., without memory) in the described ID and 2D 
settings. We compute the ratio of (t(x)) for spiders with 
r = 1 and (r(x)) for spiders with r = r pt(x). This ratio 
is plotted in Fig. [14] against the surface size for ID and 
2D. 

The plot shows that the cleavage rate r gives spiders 
more advantage on a 2D plane searching for a circle than 
on a ID strip searching for its ends. This advantage can 
be attributed to the amount of substrates that spiders 



FIG. 14: Ratio between {r(x)) of the spiders with r = 1.0 
and r = r opt . 



leave behind as they move away from the origin. In ID, 
spiders do not leave any substrates behind, so there are 
no substrates between the left and the right ends of the 
sea of products, and when a spider moves back towards 
the origin there are no substrates to bias it towards the 
boundary. In 2D, the shape of the product sea can be 
very complicated and there can be many substrates left 
behind. When a spider moves backward it has a high 
probability to still encounter substrates, which can bias 
it towards the boundary. The higher r op t(x) values in 
2D can be attributed to the direction of the emergent 
bias when the spider is on the border between visited 
and unvisited sites. In ID, the border is simple, and its 
shape remains the same over time. It is defined by the 
two closest unvisited sites to the origin on the right and 
on the left side. As a result, when the spider with rate 
r < 1 is on the border, it is always biased in the desired 
direction — away from the origin. In 2D, the shape of 
the border can be much more complex and is even not 
necessarily connected. This type of border leads to a 
much weaker bias towards the edge of the surface. In 
many cases the spider is not biased directly to the edge 
in the direction of the shortest path, and sometimes the 
spider even is biased back towards the circle's center. 
Despite that, the greater amount of substrates accessible 
by the spider in 2D overcomes the weaker bias and makes 
spiders more efficient at finding the absorbing boundary. 
Fig. US] shows the average density of products; the spider 
has a high probability of encountering substrates when it 
turns back towards the origin. 



VII. DISCUSSION 

Our simulations show that the MFPT of two-legged 
spiders depends strongly on the kinetic parameter r in 
all studied cases. The MFPT is lower for r < 1 (i.e., 
with memory) than for r = 1 (i.e., without memory) in 
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FIG. 15: Average density of products when spider with 
0.25 and k = 2 is at distance 1000. 



the one-dimensional case when the spider is searching for 
the ends of a one-dimensional track. For one-legged spi- 
ders the parameter r does not affect leading asymptotic 
behavior, however it slows them down by increasing the 
constant of the sub- leading term. In the extension of this 
problem to two dimensions, when the spider is searching 
for the contour of a circle from its center, the advantage 
of having r < 1 is even more significant, despite the less 
effective bias provided by the shape of the leftover sub- 
strates. Here the bias provided by the substrates does not 
always direct the spider towards the absorbing boundary. 
In contrast, on a ID track, the substrates always bias the 
spider towards the closest end when one of the legs is at- 
tached to them. The disadvantage in 2D is overcome by 
the greater amount of substrates that are accessible to 
spiders. In ID, the spider (of the studied gait) does not 
leave any substrates behind when it progresses towards 
the ends of the track. In 2D, the shape of the sea of 
products is complex, and many substrates are left be- 
hind. Those substrates can bias the spider towards the 
absorbing boundary when it starts to turn back towards 
the origin. 

When we reverse boundaries in 2D we make the ab- 
sorbing boundary the single site in the center of the circle, 
and start the spider from the contour, the parameter r 
can still improve the MFPT. It is difficult to find the sin- 
gle target site in this scenario, and the MFPT increases 
significantly for all types of walkers. However, it is plau- 
sible that when most of the surface is explored, spiders 
with r < 1 are more likely to stick to small remaining is- 



lands of substrates (i.e., unvisited sites). The presence of 
these islands increases the probability that spiders with 
r < 1 will find the target site, as opposed to a spider 
with r = 1, which would not react to the presence of 
those islands. This is a likely explanation for the greater 
importance of r for this boundary condition. 

Despite the varying importance of the parameter r in 
the studied scenarios, in all cases there exists an opti- 
mal value of r that minimizes the first passage time to 
the absorbing boundary. To the extent that catalysis is 
accessible as a design parameter of molecular spider as- 
semblies, the results provide a way to optimize spider 
system performance in various target search scenarios. 



Acknowledgments 

The authors would like to thank Paul Krapivsky for 
detailed discussions concerning their model and analy- 
sis. This material is based upon work supported by the 
National Science Foundation under grants 0829896 and 
1028238. We would also like to thank NVIDIA Corpora- 
tion for a hardware gift that made possible some of the 
simulations. 



Appendix A: Optimum value of r 

First, we find the derivative of the assumed form pf 
(t(x)) (Eqs.lBland lAljl with respect to r. 

d(r(x)) _ ci(c 3 r + c 4 ) - c^^r + c 2 ) 2 



dr 



(C3f + c 4 ) 2 
c 5 (c 7 r + c 8 ) - c 7 (c b r + c 6 ) 
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(c 7 r ■ 

C3C2 2 , 
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The fitting of Eq.[6]to the estimates of Ai and a\ yields 
estimates of constants c\ to eg. Next, we substitute the 
constants into Eq. IA1I and find when the derivative is 
zero. 
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